Theory of the Kagome Lattice Ising Antiferromagnet in Weak Transverse Field 
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We study the quantum Ising antiferromagnet on the Kagome lattice, with weak transverse field 
dynamics and other local perturbations. We analytically demonstrate the possibility of a disordered 
zero-temperature phase that is smoothly connected to the phase at strong transverse fields. This 
is done by means of an appropriate mapping to a compact U(l) gauge theory on the honeycomb 
lattice that is coupled to a charge-1 matter field. Our results are consistent with existing Monte- 
Carlo calculations. The differences with other commonly studied lattices (in which such disordered 
phases do not obtain at weak transverse fields) is explained. Ordered phases are also shown to be 
possible in principle in the weak transverse field limit, and are briefly studied. 

PACS numbers: 



I. INTRODUCTION 

Geometrically frustrated quantum magnets hold a very 
important role in the active and challenging search for 
exotic quantum phases and spin liquids. It is hoped 
that the strongly competing exchange interactions be- 
tween spins may in some cases prevent not only devel- 
opment of magnetic order, but any other long range or- 
der as well. In certain circumstances the resulting state 
(known as a "spin liquid" ) , disordered by quantum fluc- 
tuations, could posses various unusual properties, such as 
non-trivial topological structure and fractionalized exci- 
tations. Much of the effort to find such fractionalized 
spin liquid states has been motivated by proposals of 
their relationship to the pseudo-gap phase of high tem- 
perature superconductors,^"'^ as well as many other un- 
usual systems. ^"^. In addition, their topological proper- 
ties seem very promising for future applications in quan- 
tum computing. ^'^ 

Out of many lattices that can host a frustrated mag- 
net, the pyrochlore and Kagome attract most attention. 
Quantum liquid arising from only the nearest neighbor 
interactions has not yet been ruled out in these systems, 
unlike in many others which initially had been looked 
upon with great hope. The reason for this is most prob- 
ably their corner-sharing structure that yields an ex- 
tremely large classical degeneracy. The two-dimensional 
Kagome lattice (Fig. 1) has been a promising host for ex- 
otic physics for more than ten years. Various physical 
states have been found at low temperatures in experi- 
ments on non- ideal Kagome magnets with spin S ^ 3/2, 
ranging from magnetically ordered ones in Jarosites,^^ 
to the spin-glass in SrCrgpGai2-gpOig (SCGO)^^'^^ and 
newer Ba2Sn2ZnGazCr-j022 (QS ferrite).^"' Numerical 
computations on small samples have provided hints of 
a liquid caused by quantum fluctuations in the ideal 
Kagome lattice spin 5—1/2 Hcisenberg model. ^^"^^ Fur- 
thermore, in comparison to other spin systems, ^^'^^ this 
one has a rather unusual spectrum, consisting of gapped 
spin-carrying excitations, and a band of seemingly gap- 
less singlet states below the spin-gap. There have been 
several theoretical attempts to understand the unusual 



spectrum of the Kagome lattice quantum Heisenberg an- 
tiferromagnet, and the nature of its ground-state. Some 
of them favored a spin liquid, ^^'^°"^^ while some other 
opened up a possibility of a valence bond crystal. ^^'^^ 
Apparently, the Kagome lattice antiferromagnets are ex- 
cellent and promising systems in which exotic phases of 
matter could be found. 

A simpler theoretical context in which the effects of 
quantum dynamics on classical geometrically frustrated 
magnets may be studied is provided by considering frus- 
trated Ising models perturbed by a transverse magnetic 
field. If the transverse field is strong (much larger than 
the typical exchange energy scales of the Ising interac- 
tion), disordered paramagnetic phases will result. The 
interesting and non-trivial questions therefore arise in the 
limit of weak transverse fields. What is the effect of such 
a weak transverse field on the classical Ising magnet? 
Generally a number of different outcomes are possible. 
Consider the rather common situation where the clas- 
sical Ising magnet has a macroscopic number of degen- 
erate ground states, and no long range magnetic order 
even at zero temperature. In some such cases a weak 
transverse field is known to immediately select a partic- 
ular ordered ground state from this degenerate manifold. 




FIG. 1: Kagome lattice is a corner-sharing two dimensional 
lattice. The frustrated units are triangular plaquettes, and 
they are only minimally connected into the lattice. 
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This phenomenon, known as "ordcr-by-disorder" , hap- 
pens for instance for the triangular Ising antifcrromag- 
net, or the fully frustrated square lattice Ising model in 
weak transverse fields. In contrast, on the Kagome lat- 
tice Monte-Carlo calculations^®'^'' report that the clas- 
sical Ising antiferromagnet stays disordered even in the 
presence of weak transverse fields. This phenomenon has 
been dubbed "disorder-by-disorder". 

Our interest in this paper is in elucidating better the 
physics of the Kagome Ising antiferromagnet in weak 
transverse fields, and possibly other local perturbations 
consistent with the Ising symmetry. Our particular fo- 
cus will be on obtaining some analytical understanding 
of the disordered paramagnetic phase found in numeri- 
cal calculations. To that end we develop a reformulation 
of the model in terms of an appropriate gauge theory 
that lives on the honeycomb lattice. General structure of 
the phases of this gauge theory will be analyzed through 
duality transformations. Our approach enables us to un- 
derstand differences between the Kagome lattice and the 
analogous problems on the square or triangular lattices 
where order-by-disorder seems to occur. 

Besides its intrinsic interest, the problem studied in 
this paper also serves as a useful technical playground 
for learning about some of the fundamental issues that 
arise in analytical studies of frustrated magnets. The- 
oretical demonstrations of quantum liquid-like behavior 
in specific microscopic models of spin or boson systems 
are rare, and currently involve interactions beyond the 
nearest neighbor. The ideas developed in this pa- 
per may perhaps motivate construction of simpler mi- 
croscopic models that display exotic phenomena such as 
fractionalization. We emphasize however that the para- 
magnetic "liquid" phase discussed in the present paper 
is non-exotic, and is expected to be smoothly connected 
to that which obtains in the large transverse field limit. 

II. LATTICE FIELD THEORY OF 
FLUCTUATING SPINS 

The basic model that we want to study is the antifer- 
romagnetic Ising model in transverse magnetic field on 
the Kagome lattice (TFIM): 

H = jY,S:S^-TY,S! ■ (1) 

Here i,j label the sites of the Kagome lattice, and Si is a 
spin S = ^ moment at site i. J is the antiferromagnetic 
Ising interaction strength, and F is the strength of the 
transverse magnetic field. More generally, we will focus 
on the regime where the energy scale J, which fixes an 
easy axis, is much larger than all other energy scales (J » 
r), and consider various ways in which dynamics can 
be given to the spins, without conserving any quantities 
(clearly, the alternate limits of large T is trivial). The 
main goal is to study the structure of possible phases that 
can emerge when the frustrated Ising antiferromagnet is 



endowed with weak quantum dynamics (that preserves 
Ising, but no other spin symmetries). 

Our strategy is as follows. We derive a low energy ef- 
fective theory that is appropriate in the easy axis limit. 
This may be represented as a compact U(l) gauge the- 
ory on the honeycomb lattice that is coupled to a bosonic 
"matter" field (with gauge charge 1). There is in addi- 
tion a non-zero static background charge at each site. 
The utility of a compact U(l) gauge theory (with ap- 
propriate background charges) to describing the low en- 
ergy physics of frustrated easy-axis magnets has been 
pointed out several times in the literature. However, in 
contrast to the Kagome lattice, on other lattices typically 
the gauge theory has no dynamical matter fields. In two 
spatial dimensions on these other lattices the gauge the- 
ory is in a confined phase, and the presence of background 
charges leads to broken translation symmetry. Presence 
of the additional dynamic matter field distinguishes the 
Kagome lattice from these other lattices. As we will see, 
it is now possible to have a phase that is also "confining" 
(more precisely a Higgs phase), which preserves transla- 
tion symmetry even in the presence of the background 
charges. For the original Kagome TFIM this describes 
a translation invariant paramagnet. The general possi- 
bility of such translation invariant Higgs phases in such 
gauge theories has been discussed before. ^'^ 

We analyze the gauge theory appropriate for the 
Kagome TFIM using duality transformations. The dual 
theory will turn out to be equivalent to a certain XY 
model with a three-fold anisotropy. The disordered phase 
of this model is the Higgs phase discussed above. The 
ordered phase describes a situation where there is order- 
by-disorder in the original Kagome magnet. 



A. U(l) Gauge Theory 

The nearest neighbor Ising coupling on the Kagome 
lattice can be conveniently written as a sum of terms 
defined on the triangular plaquettes: 

{ij) A iGA 

This allows us to easily describe the sector of low energy 
states: total spin on every triangle should be ±1/2. Spin 
configurations that satisfy this condition are least frus- 
trated. Let us express these states using a set of variables 
defined on the honeycomb lattice, whose sites we will la- 
bel by p and q. Figure 2 illustrates relationship between 
the Kagome and honeycomb lattices. The honeycomb 
bonds contain Kagome sites, and we can associate with 
them the Kagome spins: = S^^^^^ . On the other hand, 
the honeycomb sites reside inside the Kagome triangular 
plaquettes. It is useful to keep track of the total spin 
(whether it is +1/2 or —1/2) on any triangular Kagome 
plaquette. We therefore introduce a variable that mea- 
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sures the total plaquette spin: 
(Vp) s 



(3) 



We impose restriction to the low energy sector by re- 
quiring that Sp take only values ±i. Note that when a 
Kagome spin is flipped (provided that flipping does not 
introduce more frustration), on the triangles that con- 
tain it change sign in the same direction. With only the 
Ising interaction present there is a large number of classi- 
cal ground states: every spin configuration that satisfies 
Equation (3) with = ±1/2 is a classical ground state. 
Inclusion of the T or other terms in the Hamiltonian will 
split this huge degeneracy of the ground state manifold. 

It is possible to perturbatively construct an effective 
theory that describes the low energy dynamics in the 
ground state manifold, and express it on the honeycomb 
lattice: 



(pg) 



I xs't + h.c. 

ipq) q 



(4) 



* ('^(12) '^{23) '^(34) '^(45) "^(56) '^{61> 
O 

+ /i.e.) . 

The Hilbert space of this theory is defined only by the 
least frustrated states. The lowest order dynamical term 
is a single Kagome spin flip caused by the transverse 
field, and the operators simply project-out the states 
that are not minimally frustrated. At higher orders of 
perturbation theory multiple spins are being flipped, and 
the smallest "ring-exchange" term appears at the sixth 
order, so that t ^ -p. Alternately, we could have imag- 
ined adding such a ring-exchange term to the original 
model (in addition to the transverse field term) . We will 
examine the properties of this effective Hamiltonian in 
the low energy subspace for arbitrary F and t. The pure 
transverse field model then corresponds to the particular 
limit r ^ All fluctuations are constrained by (3). The 
effective theory will have this general form for all kinds 
of Kagome spin models with an easy axis, provided that 
dynamics preserves only the Ising symmetry. 

The Hamiltonian (4) can be interpreted as a model 
of charged bosons moving in presence of a fluctuating 
electromagnetic field. In order to formulate this inter- 
pretation, we need to exploit the bipartite nature of the 
honeycomb lattice. Let us introduce a fixed field Sp that 
takes values -1-1 and —1 on two different sublattices, as 
shown in Fig. 3. Then, define: 



En 



(5) 
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Note that E„ 



-Ear), SO that it can be viewed as a 
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FIG. 2: Relationship between the honeycomb and Kagome 
lattices. Every honeycomb bond contains one Kagome site, 
while every Kagome triangle contains a honeycomb site. 



will be interpreted as an integer-valued "electric" field. 
The integer Up will be interpreted as the gauge charge of 
a bosonic "matter" field that couples to the electric field. 
The constraint (3) becomes: 



(Vp) ^ Epq ^Up + e 



p 1 



(6) 



qGp 



where the sum on the left-hand side is taken over three 
honeycomb sites neighboring to p. Natural interpretation 
of this equation is Gauss' Law: divergence of the electric 
field Epq is equal to the total local charge. Note that the 
boson occupation numbers Up take values and Sp. There 
is an additional fixed background charge distribution Sp. 

At this level, the boson occupation number and electric 
field strength are constrained to only two integer values 
by (5). It is useful to soften this "hard-core" by allowing 
Up and Epq to take arbitrary integer values, but penal- 
izing fluctuations where either quantity assumes values 
different from that dictated by the hard-core condition. 
It will also be useful to introduce the corresponding con- 
jugate operators, ipp and Apq, which are angular variables 
in [0, 27r). The boson creation and annihilation operators 
Sp simply become exp(±zep(^p), and similar holds for the 
electric field: S^^^^ eyi-p{±iepApq). Now it is straight 
forward to rewrite the Hamiltonian (4) as a compact 
U(l) gauge theory (up to a constant): 

h = u,Y: {Epq - ) ' + c/2 E ("p - y ) ' 

{pq) P 

- T^COS^iPq-iPp-Apq'j (7) 




vector living on the bonds of the honeycomb lattice. This 



We have labeled by E^q^ a fixed background electric field 
that originates due to the background charge (see Fig. 3). 
The terms proportional to C/i and U2 penalize fiuctua- 
tions of the boson number and the electric field away 
from the preferred "hard-core" values. 
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FIG. 3: Fixed background charge Sp = ±1 on the honey- 
comb sites, and the electric field created by it (up to a 



multiplicative factor) . 



(0) 



£p/2 



-E, 



(0) 



On every honeycomb bond we define 
= -£,/2 (divS(°) = fep). 



We see that Apq plays role of a vector potential. The 
term at the lowest order of perturbation theory describes 
boson hopping on the honeycomb lattice, while the term 
at the sixth order gives energy cost to the "magnetic flux" 
(away from 2-jt x integer values). It is also easy to write 
down other terms consistent with the symmetries and the 
U{1) gauge structure. Fluctuations arc subject to the 
Gauss' Law (6). Note that the charged bosons cannot 
screen out the background charge (rip — » {Q,£p}) with- 
out paying a large price {1/1,1/2) ■ This is a consequence of 
magnetic frustration in the original Kagome spin model. 
In fact, there is a large number of degenerate boson con- 
figurations that accomplish the best possible screening, 
and they correspond to the least frustrated states. How- 
ever, the very fact that there is a matter field in this com- 
pact U(l) gauge theory distinguishes the Kagome lattice 
from other commonly studied lattices. It gives hope that 
a translation symmetric paramagnetic phase could exist 
on the Kagome lattice, in contrast to the situation for the 
triangular Ising antiferromagnet, or the fully frustrated 
square lattice Ising magnets (with weak transverse field 
dynamics). 



B. Duality transformation 

In this section we perform a duality transformation on 
the compact U{1) gauge theory derived above. This will 
enable us to analyze the structure of the possible phases. 
The duality transformation proceeds in a standard fash- 
ion. We first derive the path-integral form of the compact 
U(l) gauge theory. All terms denoted by ellipses in (7) 
will be ignored. The action will contain a usual Berry's 
phase (we will omit the time index): 



r \{pq) p J 



(8) 



and a potential energy part: 



U,J2{Epq-Ell^)\u2j:{np-'-^y 
(pq) P 



where 5t is the imaginary time increment. The kinetic 
energy part will be obtained by applying Villain's ap- 
proximation to the cosine terms in (7): 



cos 9 



E 



m— — 00 



-Km —imO 



t = 2e 



-K 



0. (9) 



Two new integer-valued Villain fields will appear: a par- 
ticle current jpg, and a magnetic field scalar Br that lives 
inside plaquettes of the honeycomb lattice, or equiva- 
lently on the dual triangular lattice sites (see Fig. 4): 



2 

pq 



(10) 



Or 



{pq) 



+ i^jpq {V>q -^p- -Apq) +i^Bri^Apq 
{pq) Or \{pq) J 

We will treat the constants K\ and as free param- 
eters (they are in principle determined in terms of the 
microscopic parameters that define the original Kagome 
Hamiltonian) . Our interest is in exploring the general 
nature of the possible phases that are contained in this 
dual action. The angular variables ^pp and Apq can now 
be formally integrated out, yielding Kronecker-delta fac- 
tors in the path-integral for all integer- valued expressions 
that they couple to. Such factors simply express famil- 
iar laws of electrodynamics. The integral over the boson 
phase angle will give rise to the current conservation law: 



2-n- 



dtpp 



Knp^^jpq 

qep 







(11) 



while the integral over the vector potential will reproduce 
a two-dimensional Maxwell's equation: 



p2tt 
I ^Apq 

Jo 



A^-Epq + jpq = Br — Br' ■ (12) 



The direction of the triangular lattice vector Br — Br' in 
the last equation is related to the direction of Epq by the 
right-hand rule (see Fig. 4). 

The equations (6), (11), and (12) can be solved on the 
dual triangular lattice, whose sites will be labeled by r 
and r'. As is standard, we first define spatial components 
of a dual gauge field Arr' on the dual triangular lattice, 
and a fixed background gauge field A^^J, , such that their 
circulations on a triangular plaquette are determined by 
the charge contained "inside" the plaquette: 



(Vp) 



E- 

{rr') 



E 

{rr') 



A 



(0) 



(13) 
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that remained after integrating out the conjugate angles: 
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FIG. 4: Duality between the honeycomb and triangular lat- 
tices. The triangular lattice sites sit inside the honeycomb 
plaquettes, and vice versa. Duality between a honeycomb 
lattice vector p ^ q and a triangular lattice vector r ^ r' 
is shown in the lower-left portion of the graph: their direc- 
tions are related by the right-hand rule. Divergence on the 
honeycomb lattice translates into negative lattice curl on the 
triangular lattice (we always take circulation in the counter- 
clockwise sense). 



The background gauge field A'f^, can be specified in many 
different ways, and we will discuss a convenient choice 
later. For now, we will only assume that it does not 
depend on imaginary time. If we substitute these defini- 
tions into the Gauss' Law (6), and note that the hon- 
eycomb lattice divergence of the electric field translates 
into the negative triangular lattice curl, we can solve it 
by writing: 



Xr - Xr 



A 



(0) 



(14) 



Again, orientation of the vectors on the right-hand side 
is related to that of Epq by the right-hand rule. The new 
field Xr is an allowed degree of freedom, since lattice curl 
of a pure gradient vector field is zero. Care must be taken 
to ensure the integer- valued nature of Epq. A convenient 
way to enforce this will be discussed below. Let us also 
define temporal components of the dual gauge fields: 



(Vr) Ar^r+f = Br - ArXr 



A, 



r,r+r 



(15) 



Substituting this and (14) into (12) gives us an expres- 
sion for the particle current: 



Jpq 



A 



r,r-fT 



ArAr 



■ curU . (16) 



Therefore, the current becomes expressed as lattice curl 
of the dual vector potential, taken on the triangular lat- 
tice temporal plaquette that is pierced by the current 
vector j = jr. Finally, we note that the current conser- 
vation (11) is not independent from the Gauss' Law (6) 
and the Maxwell's equation (12). 

We can now obtain dual form of the following action 



(pq) 



(17) 



We will eliminate Epq, Br, rip, and jpq using (14), (15), 
(13), and (16) respectively. The particle number and 
current terms together yield curls of the dual vector po- 
tential on spatial and temporal plaquettes of the trian- 
gular lattice. Though not necessary, for simplicity we set 
U25t = K2 = g/2, and UiSt = Ki = e^/2 and consider 
the phase diagram in the resulting section of coupling 
constant space. In order to complete the duality trans- 
formation, we must also translate the quantities Ep^q\ 
and Ep into the dual language. A natural dual counter- 
part of the background electric field vector Ep^ is A^^J, , 
since both are determined by the background charge dis- 
tribution Ep. In the spirit of equation (14) this suggests 



the identification E^^ = 



(0) 



2 rr 



shown in the Fig. 5. 

However, the values of A^^J, are then not integers, and 
compensation is necessary in order to keep Epq integer- 
valued in (14). A convenient solution is to require that 
Xr fields take particular non-integer and site-dependent 
set of values. Specifically, we demand that Xr—Xr'^ be in- 



tegers, with the fixed fractional offsets x 
in the Fig. 6. Then, the dual action is 



(0) 



as illustrated 



S='^j:{xr-Xr--Arr- + lA^:jy 

{rr') 

-f|5:(curl(A-iA(°)))', (18) 



plaq. 



where the first summation extends over all space-time 
links, and the second over all spatial and temporal pla- 
quettes. Let us also shift 



A--A0 . 
2 " 



(19) 



The action then reads: 

„2 



S = 



{rr') 



Xr - Xr 



Xr J 



Arr') +1 ^(curU)' . (20) 

plaq. 

and Arr' + take inte- 



are given in the 



In this theory Xr 

ger values. The fractional residua Xr 
Fig. 6(d). 

It will actually be convenient to impose these integer 
conditions on the Xr and Arr' fields softly (in the same 
spirit as the usual sine-Gordon description of Coulomb 
gases). The result is a generalized sine-Gordon theory 
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FIG. 5: Convenient definition of tfie fixed background gauge 
field On every bond Af), = 1/3 in the given direction 

r — »■ r', so that its (counter-clockwise) circulations yield Ep on 
every plaquette. This vector field is formally dual to that of 
E^^ on the honeycomb lattice (up to a factor). 




(a) 



(b) 




(c) 



(d) 



FIG. 6: (a) Fixed background gauge field (b) A curl- less 

vector field A'^^,. Every arrow contributes 1/3 in the given 
direction, so that the sum A^^^, H- A'^^, is an integer-valued 
vector field, (c) A scalar field Xr whose gradient gives the 
curl- less vector field: A'^^, — x'r' " x'r- (d) Fractional parts 
Xr°'' of the scalar field x'r- Requiring that the fields Xr in (14) 
have fractional parts equal to Xr"^ compensates for fractional 
values of A^°\, and makes Epq an integer. 



with the structure: 

c2 



7 (xr - Xr' - Arr'^ | ^ (curU) ' 

(rr') plaq. 

i^5^cos27r(A.' + ^4°') 



(rr') 

-7^cos27r(xr-x[°^) - 



(21) 



In principle, higher harmonics of the lowest order cosine 
terms shown above must also be included. The fields Xr 
and Arr' now take real values. Let us absorb the field Xr 
into Arr' by the shift: 



+ Xr - Xr' 



(22) 



to get the action: 



{rr') plaq. 

- COs27r(xr - Xr' + Arr' + ^4r') 

(rr') 

- 7^cos27t(xr-x[°0 • 



(23) 



Note that the term proportional to appears as a 
"mass" term for the gauge field Arr' ■ We therefore inte- 
grate out Arr' ■ This may be explicitly done in the follow- 
ing manner. We formally expand exp(— 5*) in powers of 
K, and decompose every cosine factor from the expansion 
using the Euler's formula 2cos^ = exp(i6') -|- exp(— 16*). 
The expansion takes the following form: 



e-^ = exp{7^cos2^(xr~xl°0} 



X E C{r^ 
{m„,} 

, 9 



(24) 



{rr') 



plaq. 

+ E 27riTOrr' (xr " Xr' + Arr' + ^^tr') } ' 
r 

where integers, and denotes their dis- 

tribution on the whole lattice. This is a Gaussian form, 
and Arr' can be easily integrated out. Clearly, this causes 
{mjr'} dependent renormalization of the C{rn^^,} factors, 
which in turn corresponds to renormalization of K. After 
resummation over mrr' , an XY model is obtained at the 
lowest order in (renormalized) K: 



S = —K cos27r^Xr 



Xr 



{rr') 



2 



(0) 



(25) 



-7Ecos27^(x.-X^°^)-0(^^7')• 
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The variables I-kXt should be treated as angles. Adding 
integers to Xr in the gauge theory (20) can always 
be compensated by a gauge transformation, so that the 
gauge inequivalent states in (25) correspond to different 
fractional parts of Xr- This is why the obtained effective 
theory is an XY model. Physically, the dual field 27rxr 
represents the phase of an operator that creates 27r units 
of vorticity in the bosonic matter field of the original hon- 
eycomb lattice U(l) gauge theory (7), and thus must 
fundamentally be a phase, after the redundant gauge de- 
grees of freedom have been removed. The structure of 
this theory is such that the XY field fluctuates in pres- 
ence of a fixed staggered flux, and there is also an external 
field that apparently explicitly breaks the XY (and ap- 
parently also lattice) symmetries (see Fig. 7) . The physi- 
cal meaning of this explicit symmetry breaking term is as 
follows. In the original honeycomb lattice gauge theory 
27r vortices of the bosonic matter field carry 27r units of 
the U(l) gauge flux. As the theory is compact, instanton 
events where this gauge flux changes by 27r arc allowed - 
thus the vorticity is not strictly conserved. The explicit 
XY symmetry breaking term in the dual representation 
precisely describes these instanton events. 

What is the actual symmetry of this action? As dis- 
cussed above, when 7 = there is a global XY symmetry 
that is apparently broken when 7 7^ 0. However, a dis- 
crete Z3 subgroup of the global XY symmetry survives 
when combined with translation by one lattice spacing. 
For instance, the action is invariant under translation 
along the horizontal a;-axis by one unit: 

2 

Xr Xr-x + 2 • (26) 

Similar transformation properties obtain under other 
translations (as well as the other lattice symmetries). 

C. Phase Diagram 

Armed with the dual formulation, we can now see qual- 
itatively why a translation invariant phase is possible in 
this model. Consider the limit 7 = 0. Then the resulting 
global XY model will (for K small enough) possess a dis- 
ordered phase with short-ranged correlations for the Xr 
field. Upon increasing K a transition to an ordered phase 
with some pattern of XY ordering will occur. Now con- 
sider turning on a small 7. Its effects will be innocuous in 
the small-ii' disordered phase. In particular, the discrete 
Z3 symmetry of | shifts of the Xr field, which realize 
lattice translations (26), will continue to stay unbroken. 
This is therefore a translation-invariant phase. It is read- 
ily seen that this phase is also invariant under all other 
lattice symmetries. The 7 term will have a much more 
significant effect in the large- ii' ordered phase, where the 
XY symmetry is completely broken: it will pin the overall 
orientation of the XY ordering pattern. Lattice symme- 
tries are consequently expected to be broken in such a 
phase. 




(a) (b) 

FIG. 7: (a) Field lines of the background gauge field A^^^,l2 

that produces a staggered flux, (b) External XY field 27rxr°' 

that breaks XY and translational symmetries (compare with 
Fig. 6(d)). This field is constant on any given sublattice of 
the triangular lattice. 



The existence of the disordered translation-invariant 
phase is our primary conclusion. How do we think about 
it in terms of the original gauge theory? From the phys- 
ical discussion above, the field e^^^^^ creates 2tv flux of 
the original gauge field, which in turn is bound to a 27r 
vortex in the phase of the bosonic matter field. Thus 
the disordered phase is to be thought of "Higgs" 
phase where the bosonic matter fields have condensed - 
which gaps out their vortices. Indeed, our gauge model 
is closely related to a similar one discussed in Ref.^^ (the 
" A?' = 1 S J model" ) where similar phenomena were shown 
to arise. Note that the Higgs phase is preferred by the 
boson hopping term in the gauge theory. Thus it is rea- 
sonable to expect that this phase is realized in the limit 
T t. This expectation is indeed consistent with the 
numerical results of Ref.^^'^''. 

We can formally back up this discussion by considering 
an alternate soft-spin version of the dual XY model. To 
that end we define: 

$r ~ e^""'^^ . (27) 

We have labelled the triangular lattice sites by their cor- 
responding vectors r. After substituting this into the XY 
model (25) and relaxing the "hard-spin" condition on $r 
we obtain a "soft-spin" lattice theory: 

(rr') 

+r\^r\^ + U\^r\'^ + ■ ■ ■ (28) 
-7^(e-^«-$, + /i.c.)---- . 

r 

The quadratic part of this action can be diagonalized to 
obtain: 

Y,{-Ke^ + rm^W (29) 
q 
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(a) 



(b) 



FIG. 8: (a) Direct triangular lattice, and a choice of primitive 
vectors: ai — x,a2 = \x + -^y. (b) Reciprocal lattice with 
primitive vectors 61 = 2-k(x — "75^)1 ''2 = 27r ■ -^y- Shaded 
region represents the first Brillouin zone. Emphasized corners 
of the Brillouin zone represent the wavevector Q = (three 
equivalent points). 



where 



Eq = COS I 9, 



COS 



+ COS 



3/ 



(30) 



It is straight-forward to show that Sq is maximized at 
two difFerend wavevectors: q — 0, and q — —Q, where Q 
is the wavevector that describes spatial variation of the 
external XY field: 



^ 47r , 



(31) 



In order to aid understanding of this special wavevector, 
we plot its location on the reciprocal lattice in the Fig. 8. 

In the following, we focus on the limit K The 
fluctuations at low energies will be dominated by the 
modes in vicinity of the two different wavevectors, q = 
and q = ~Q. Such fluctuations can be expressed by: 



^1 



-iQr 



^^2 



(32) 



where V'l.r and V'2.r are fields that vary slowly on the 
scale of the lattice spacing. It is then useful to go to 
a continuum limit which focuses on the long wavelength 
fluctuations of these two fields. This continuum theory 
must be expressed in terms of ipi^r and "02, r only, and 
must be invariant under all symmetry transformations of 
the lattice theory (28). First, we note that the quadratic 
part (29) is invariant under global XY rotations $r 
e'^$T-, as well as all lattice symmetry transformations at 
small deviations k from the wavevectors q G {0, — Q}: 



(Va) 



const. + Kfc2+0(fc3) . (33) 



In addition, there is also a symmetry under exchanging 
the two kinds of slowly varying fields, V'l.r "02, r- On 
the other hand, as discussed in the previous section, the 



7 term (and various similar higher order terms) break the 
XY symmetry explicitly, and reduce the symmetry group 
to a discrete set of transformations. It is easy to show 
that under a unit translation by r the fields transform 
as: 



iQr 



^2 



-iQr 



1p2 



(34) 



Note that this rotates the phases of ipi and tp2 by 120 
degrees in opposite directions. Similarly, the symmetry 
transformation ipi^r ^ "02, r of the quadratic part of the 
action must be accompanied by lattice inversion r —r 
in order to keep the whole action (28) invariant: 



V'l.r i>2 



0: 



2,r 



(35) 



The remaining lattice transformations do not give rise 
to further reduction of the symmetry group. We may 
now write down a Landau-Ginzburg effective field theory, 
symmetric under transformations (34) and (35): 



dfc 



(a + bk^ 



dr 



+ h.c. 



(36) 



/3(<rV'2,r +V^2,r^l,r- + /l-C-) + 
l{Kr+^lr + h-c) 

{\A,r\'+\^2j') +u'\A,r\'\^2,r 



dr 



l" (V'l.r '/'2,r- 
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l,r I 



l^2,r 



It is clear that at least for the coupling r sufficiently large 
and positive a stable disordered phase will exist where the 
Z3 symmetry is unbroken. This then corresponds to the 
translation symmetric phase of the original lattice model. 

Ordered phases are also possible in this theory. In or- 
der to reveal their structure, we simply consider the static 
classical XY field configurations that minimize the action 
of the effective XY model (25). Let us ignore the terms 
0(iir^,7^), and restrict our attention to configurations 
where the value 27rxr depends only on the sublattice to 
which the site r belongs. The energy per unit-cell of such 
configurations is: 



E 



-3K 



cos2^(xi-X2 + ^) (37) 

-hcos27r(x2 - X3 + ^) 

-f cos27r(x3 ~ + ^) 
3 

7^cos27r(x„-xi°0 



where we have labeled the sublattices by 1,2,3 in such 
a way that the background gauge field vector A^^-*, cir- 
culates in the direction 1 ^ 2 ^ 3 ^ 1. Two kinds of 
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states can be found. For sufficiently large 7, the XY field 
simply follows the "external" XY field: Xr = X-r \ and 
an XY "spin density wave" is established at the wavevec- 
tor Q. Naively, the XY rotation and lattice symmetries 
appear explicitly broken in this state. However, the an- 
gular order parameter 2TrXr = 277X7-°' is invariant under 
the lattice translations (26) (and other lattice transfor- 
mations), so that this state of the effective XY model 
corresponds to the disordered state of the original spin 
model. On the other hand, if 7 = the ordering will be 
determined by the nearest-neighbor interaction K. As 
we have argued before, there are two ordering wavevec- 
tors {q = and q = —Q) preferred by the K term, but 
neither of them coincides with the wavevector that de- 
scribes spatial variations of the external XY field Xr ''- 
Consequently, any small non-zero 7 will introduce frus- 
tration, and deform the spontaneously ordered XY "spin 
density wave" preferred by the K term. Even though 
the precise description of the ordering pattern is compli- 
cated, it is at least apparent that the XY fields Xr align 
with the external field Xr'^ 011 one sublattice of the trian- 
gular lattice, while on the other two sublattices they cant 
toward the external field, simultaneously trying to pre- 
serve the preferred ordering wavevector. Arbitrary choice 
of sublattice for the alignment gives rise to a three-fold 
degeneracy, so that the 7 term breaks the continuous XY 
symmetry down to a Z3 subgroup, associated with lattice 
transformations (as evident in the Equation (34)). How- 
ever, choice of the "parent" ordering wavevector (q = 
or q = —Q) for sufficiently small 7 is also available, and 
we will briefly discuss its physical origin. 

One symmetry transformation that we have ignored 
so far is the global spin-flip in the original spin model on 
the Kagome lattice: Sf ^ —Sf. It is straight forward to 
trace back how this transformation affects the quantities 
of the U(l) gauge theory, and its dual theory on the trian- 
gular lattice. For example, at the level of Equation (18), 
which describes the dual theory on the triangular lattice 
with integer- valued gauge field Arr', and integer- valued 
Xr — Xr^^ ) the global spin-ffip corresponds to: 

Arr'^{A^:J,+X^:!^-X^^^)-Arr' 

Xr -Xr - Xi°^ ■ (38) 

Note that the action (18) is invariant under this transfor- 
mation, and that the integer constraints are not affected. 
Only the second part of this transformation (involving 
the Xr field) survives in the effective XY theory (25) . If 
we now turn to the continuum limit, we find from (31), 
(32) and (38) that the global spin-flip is represented by: 

Since the fields V'l.r and '02,7- describe the XY "spin den- 
sity waves" at wavcvectors q = and q = Q respec- 
tively, we see that the spin-flip formally exchanges these 
two kinds of order in the effective theory. Apparently, 
this transformation is completely independent from lat- 
tice translations. Therefore, the choice of the "parent" 



r 



D 




FIG. 9: Schematic phase diagram of the effective XY model 
(25). D corresponds to the disordered phase of the original 
Kagome spin model. LRO is a spontaneously long-range or- 
dered and magnetized phase. 

ordering wavevector (q = or q = Q) in the sponta- 
neously ordered phase must correspond to the choice of 
direction of the global magnetization. 

We can finally sketch the phase diagram of the ef- 
fective XY model (25), as shown in the Fig. 9. The 
phase D is invariant under lattice translations (26) and 
global spin-flip (38), so that it corresponds to the dis- 
ordered phase of the original spin model (4), realized 
for T ^ t. The spontaneously long-range ordered phase 
LRO is obtained in the opposite limit of strong hexagon 
ring-exchange T <^ t. We argued that this phase has 
global Ising magnetization, while the translation symme- 
try is broken in a three-fold degenerate manner. This 
naturally suggests a microscopic description of the order 
parameter given in the Fig. 10. We would like to note 
that such a state exactly corresponds to the state with 
one-third of saturated magnetization on the Kagome lat- 
tice that was found responsible for plateaus in the mag- 
netization curves of some Kagome-based systems.^® The 
phase transition is most likely of the second order (if 
we had considered a small longitudinal field in the orig- 
inal Kagome lattice spin model, the ultimate continuum 
limit would have been an XY model with a three- fold XY 
anisotropy, which has a second order transition). 

III. DISCUSSION 

The theory of the Kagome lattice quantum Ising model 
that we have presented in this paper reveals that a dis- 
ordered ground state is found for weak transverse field. 
This disordered phase breaks no symmetries and is not 
topologically ordered either. It is therefore expected to 
be smoothly connected to the completely imcorrelated 
phase at large transverse field. The same conclusion has 
been strongly suggested by Monte Carlo simulations, ^^'^^ 
where no phase transition was detected as the strength 
of transverse field was varied. 
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FIG. 10: The most symmetric pattern of three-fold 
translational symmetry breaking on the Kagome lattice. 
Roughly speaking, the six spins on the emphasized hexagons 
are alternating, and coherently resonating as a singlet 
IT-lT-lT-lHHiTiTiT) ^ ^jiijg ^j^g remaining spins are ferromagneti- 
cally aligned with each another and break the global spin-flip 
symmetry. Note that in the Hamiltonian (4) the energy is 
reduced by 0{t) on every resonating hexagon. 

Such disordered phases seem to be exceptions rather 
than the rule for a large set of two-dimensional frus- 
trated lattice models with similar structure. If the clas- 
sical ground state manifold is macroscopically degener- 
ate in a system with discrete degrees of freedom, a typi- 
cal situation is that quantum fluctuations reduce the de- 
generacy down to the one associated with broken lattice 
symmetries, creating "ordcr-by-disorder" . Examples in- 
clude the hard-core quantum dimer and Ising models on 
most studied simple lattices,^^'^''. The analysis in this 
paper provides a route to understanding the differences 
between these common situations with order-by-disorder 
and the Kagome quantum Ising antiferromagnet. We 
have demonstrated how the Kagome lattice quantum 
Ising model can be described by a compact U(l) gauge 
theory with some fixed background charge at each lattice 
site. Such a gauge theoretic description has been useful 
in studies of similar Ising models on other frustrated lat- 
tices. The crucial distinction between the Kagome and 
other common systems is in the fact that the U(l) gauge 
theory of the Kagome system contains a dynamical mat- 
ter field. Without a matter field, these gauge theories 
of the frustrated magnets in two-dimensions always ul- 
timately live in a confined phase that breaks translation 



symmetry. The latter is caused by the fixed background 
charge in the gauge theory. However, in the presence of 
dynamical matter fields a translation invariant ( "Higgs" ) 
phase is generally possible. In situations where the mat- 
ter field has gauge charge 2 (as happens in the gauge 
theoretic description of quantum dimer models on non- 
bipartite lattices) such a Higgs phase also possesses topo- 
logical order and associated "vison" excitations: for in- 
stance, there is a non-trivial ground state degeneracy on 
topologically no n- trivial manifolds. In the problem dis- 
cussed in this paper the matter field had gauge charge 1. 
The resulting Higgs phase, while translation invariant, 
is topologically trivial. This then lends strong support 
to the conjecture in Ref.^^ that weak transverse fields 
on the Kagome Ising magnet lead to a phase that is 
smoothly connected to the trivial paramagnet that ob- 
tains at strong transverse fields. 

In fact, topological triviality of the discovered disor- 
dered phase is a very interesting detail. Prom the begin- 
ning of the quest for interesting quantum spin liquids, 
the frustrated magnets had been looked upon with great 
hope. In particular, the systems with extremely large 
classical ground state degeneracy have attracted consid- 
erable attention. The pyrochlore and Kagome lattices are 
thought to be ideal candidates for the spin liquid, because 
their corner-sharing structure provides such extremely 
large degeneracy, which is then dramatically lifted even 
by weak quantum fluctuations. However, in the case of 
the Kagome lattice transverse-fleld quantum Ising model, 
the corner-sharing geometry seemingly makes the lattice 
virtually disconnected. Even though there is no long- 
range order in the ground state, the obtained disordered 
phase is not a topologically ordered spin liquid. One may 
then speculate that a frustrated, but somewhat more con- 
nected lattice is probably a better platform to seek topo- 
logically non-trivial spin liquids. Another question is if a 
more correlated spin dynamics, such as the one found in 
anisotropic easy-axis Heisenberg models, is more likely to 
yield such interesting spin liquids. These issues are left 
open for future work. 
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